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SPHERICAL CONFORMAL PARAMETERIZATION OF GENUS-0 
POINT CLOUDS FOR MESHING 


GARY PUI-TUNG GHOI, KIN TAT HO AND LOK MING LUI 

Abstract. Point cloud is the most fundamental representation of 3D geometric objects. Ana¬ 
lyzing and processing point cloud surfaces is important in computer graphics and computer vision. 
However, most of the existing algorithms for surface analysis require connectivity information. There¬ 
fore, it is desirable to develop a mesh structure on point clouds. This task can be simplified with the 
aid of a parameterization. In particular, conformal parameterizations are advantageous in preserving 
the geometric information of the point cloud data. In this paper, we extend a state-of-the-art spheri¬ 
cal conformal parameterization algorithm for genus-0 closed meshes to the case of point clouds, using 
an improved approximation of the Laplace-Beltrami operator on data points. Then, we propose an 
iterative scheme called the North-South reiteration for achieving a spherical conformal parameteriza¬ 
tion. A balancing scheme is introduced to enhance the distribution of the spherical parameterization. 
High quality triangulations and quadrangulations can then be built on the point clouds with the aid 
of the parameterizations. Also, the meshes generated are guaranteed to be genus-0 closed meshes. 
Moreover, using our proposed spherical conformal parameterization, multilevel representations of 
point clouds can be easily constructed. Experimental results demonstrate the effectiveness of our 
proposed framework. 

Key words. Mesh generation. Triangulation, Quadrangulation, Spherical conformal parame¬ 
terization, Surface reconstruction. Point cloud. Multilevel representation 


1. Introduction. Contemporary scanning technologies enable efficient acquisi¬ 
tions of 3D objects. Using modern 3D scanners, data points are sampled from the 
surfaces of 3D objects for further analyses and usages. Point clouds are widely ap¬ 
plied in computer graphics, vision and many other engineering fields. However, the 
data points acquired by laser scanners are often complex and unorganized. Moreover, 
the absence of the connectivity information in point cloud data poses difficulties in 
understanding the underlying geometry of the 3D objects. This largely hinders the 
applications of the data. For instance, many applications in 3D printing [31117] and 
texture mapping [39l[24] are built upon mesh structures. With the rapid development 
of the computer industry, finding a high quality meshing framework for point cloud 
data is increasingly important. 

One possible approach for mesh generation on point clouds is to parameterize a 
point cloud to a simpler domain with the corresponding genus, such as the unit sphere 
for genus-0 point clouds. Then, a triangulation or a quadrangulation can be created 
on the parameter domain instead of the original complicated point cloud. Finally, 
a mesh structure on the point cloud can be defined with respect to the structure 
on the parameter domain. The major difficulty of computing parameterizations of 
point-set surfaces is the extremely limited information they can provide. Most of the 
existing surface parameterization methods are developed on meshes only. In other 
words, besides the locations of the point data, a given connectivity is also required 
as an input. The connectivity information plays an important role in representing 
the surface structure as well as in approximating continuous operators to minimize 
certain distortions. As a result, most conventional mesh parameterization approaches 
fail to work on point clouds. Without the connectivity information, the underlying 
geometry of the point cloud data become more obscure. Hence, it is more challenging 
in developing parameterization schemes with good quality for point cloud data. 

A good parameterization scheme of point cloud must satisfy certain criteria. In 
particular, it should retain the geometric information of a point cloud as complete as 
possible. In our case, one of the ultimate goals is to create a triangulation for a point 
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cloud by finding a Delaunay triangulation on a simpler parameter domain. It is note¬ 
worthy that in general, a mesh structure with good quality on the parameter domain 
does not necessarily imply that the associated mesh structure of the original data 
points is satisfactory. In other words, meshing the parameter domain may provide 
meaningless results if the parameterization scheme is arbitrarily chosen. Note that 
the regularity of the mesh structures is related to the angle structure of the triangles 
and quadrilaterals. To ensure the regularity of the associated mesh structure on the 
point cloud, the parameterization should preserve the angle structure of the triangles 
and quadrilaterals on the parameter domain. This motivates us the use of conformal 
mappings. 

For smooth surfaces, it is well known that the conformal parameterizations pre¬ 
serve angles and hence the local geometry of the surfaces. It is natural to consider 
the discrete analog of conformal parameterization for point cloud data. Since data 
points are sampled from real 3D surfaces, we can assume that every point cloud has 
an underlying geometry. Based on this important assumption, we consider finding 
conformal parameterizations of genus-0 point clouds. In [6], Choi et al proposed a 
fast spherical conformal parameterization algorithm for genus-0 closed surfaces in two 
steps. In the first step, a Laplace equation is solved on a planar triangular domain 
and the inverse stereographic projection is applied to obtain an initial spherical pa¬ 
rameterization. In the second step, quasi-conformal theories are applied to enhance 
the conformality of the spherical parameterization. The computation is linear and the 
conformality distortion of the parameterization is minimal. However, the algorithm 
is developed on triangular meshes only. In this work, we extend and improve the 
algorithm for point clouds with spherical topology. 

The aforementioned algorithm in [6] developed on meshes involves solving a 
Laplace equation. To extend the algorithm for point clouds, we propose a new weight 
function for enhancing the accuracy of the approximation of the Laplace-Beltrami 
(LB) operators on point clouds. Using our improved approximation, the LB operator 
in the mentioned algorithm can be accurately computed on point clouds. Also, we re¬ 
place a key step of solving for a quasi-conformal map in the mentioned algorithm by an 
iterative scheme, called the North-South reiteration^ for improving the conformality of 
the parameterizations. Furthermore, we introduce a balancing scheme for enhancing 
the distribution of the parameterization results. Experimental results demonstrate 
the effectiveness of our proposed parameterization algorithm for genus-0 point clouds. 
Our algorithm achieves global spherical parameterizations with minimal conformality 
distortions. Furthermore, with the aid of our parameterization schemes, we can eas¬ 
ily generate high quality triangulations and quadrangulations on point clouds. The 
meshes generated are guaranteed to be genus-0 closed meshes. Moreover, multilevel 
representations of the point clouds can also be easily computed with the aid of our 
spherical parameterization scheme. 

The rest of the paper is organized as follows. The contribution of our work is 
highlighted in Section In Section we review the related previous works on point 
cloud parameterizations and approximations of differential operators on point clouds. 
In Section we introduce the mathematical background of our work. In Section 
we review a spherical conformal parameterization scheme for triangular meshes, 
which is closely related to our proposed framework for point clouds. In Section we 
explain our proposed framework for spherical conformal parameterization and mesh 
generation of point clouds. In Section we demonstrate the effectiveness of our 
proposed framework by numerous experiments. The paper is concluded in Section 
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Methods 

Topology 

Parameter 

domain 

Local/global 

parameterization? 

Distortion to 
be minimized 

Meshless 

parameterization pdH [TT] 

Disk 

topology 

Plane 

Global 

/ 

Meshless 

parameterization for 
Spherical Topology [18] 

Genus-0 

Planes 

Local 

/ 

Spherical embedding [50] 

Genus-0 

Sphere 

Global 

Stretch 

Discrete one-forms [44] 

Genus-1 

Planes 

Local 

/ 

As-rigid-as-possible 

meshless 

parameterization [48] 

Disk 

topology 

Plane 

Global 

ARAP 

Meshless quadrangulation 
by global 

parameterization [62] 

Unre¬ 

stricted 

Plane 

Global 

Gradient and 
principal 
fields 


Table 2.1 

Several previous works on meshing point clouds using parameterization. 


2. Contribution. In this work, we propose a framework for meshing using 

spherical conformal parameterizations of genus-0 point clouds. Our proposed method 

is advantageous in the following aspects: 

(i) We extend and improve the spherical conformal parameterization algorithm on 
meshes in [6] for point clouds. An accurate approximation of the Laplace- 
Beltrami operator is achieved using the moving least square method [63l [25] |26] 
together with a new Gaussian-type weight function. A key step of the param¬ 
eterization algorithm in [6] for computing a quasi-conformal map is replaced 
by solving a Laplace equation on the complex plane, followed by an iterative 
scheme called the North-South (N-S) reiteration. Also, the point distribution of 
the parameterization is enhanced by a balancing scheme for point clouds. 

(ii) Our spherical parameterization method is efficient and robust to complex geo¬ 
metric structures. The algorithm completes within a few minutes and can handle 
highly convoluted point clouds. 

(hi) Unlike most of the existing approaches, our algorithm specifically minimizes the 
conformality distortion of the parameterizations. Since the local geometry is pre¬ 
served under the global spherical conformal parameterizations, we can create an 
almost-Delaunay triangulation on a point cloud by computing a Delaunay trian¬ 
gulation of its spherical conformal parameterization. The resulting triangulation 
on the point cloud preserves the regularity of that on the parameterization. 

(iv) High quality quad meshes can also be generated using our spherical conformal 
parameterization scheme. 

(v) Unlike the conventional approaches for meshing, our method is topology pre¬ 
serving. The meshes produced using our proposed framework are guaranteed to 
be genus-0 closed meshes. No post-processing is required. 

(vi) Our method is stable under geometrical and topological noises on the input point 
clouds. 

(vii) With the aid of our spherical conformal parameterization scheme, multilevel 
representations of genus-0 point clouds can be easily constructed. 
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3. Previous Works. In this section, we describe some previous works closely 
related to our work. 

Surface parameterization has been extensively studied by different research groups. 
For surveys on surface parameterization methods, please refer to mnans[42]. in 
particular, conformal parameterizations have been well established on meshes. For 
the recent works, Lai et al [21] proposed a folding-free spherical conformal mapping 
scheme by the harmonic energy minimization. In [56], Aflalo et al. proved theo¬ 
retical bounds of the conformal factor and proposed a method that minimizes the 
area distortion and avoids numerical errors of the conformal mapping. In [S, Choi et 
al developed a linear algorithm for spherical conformal parameterizations of genus-0 
closed meshes. 

In the last few decades, numerous studies have been devoted to the parameter¬ 
ization of point cloud data. In unmi], Floater and Reimers proposed the meshless 
parameterization method for unorganized point sets. The point sets are parameter¬ 
ized onto a planar domain by solving a sparse linear system. In [50], Zwicker and 
Gotsmann presented a parameterization approach for a genus-0 point cloud using a 
/c-nearest neighborhood graph of the point cloud, followed by a spherical embedding 
method for planar graphs. In ID HIE], Azariadis and Sapidis introduced the notion 
of dynamic base surfaces and suggested a parameterization scheme by orthogonally 
projecting a point cloud onto the dynamic base surface. Guo et al. m computed 
a global conformal parameterization of point-set surfaces, based on Riemann sur¬ 
face theory and Hodge theory. In [44], Tewari et al. proposed a doubly-periodic 
global parameterization of point cloud sampled from a closed surface of genus 1 to 
the plane, with the aid of discrete harmonic one-forms. Wang et al. [45] suggested 
a parameterization method for genus-0 cloud data. A point cloud is first mapped 
onto its circumscribed sphere, then the sphere is mapped onto an octahedron and 
finally unfolded to a 2D image. In [48], Zhang et al. presented an as-rigid-as-possible 
parameterization approach for point cloud data. A point cloud with disk topology is 
mapped onto the plane by a local flattening step and a rigid alignment. In [25], Liang 
et al. constructed spherical conformal mappings of genus-0 point clouds by adapting 
the harmonic energy minimization algorithm in [21] . Meng et al. m proposed a neu¬ 
ral network based method for point cloud parameterization. An adaptive sequential 
learning algorithm is applied to dynamically adjust the parameterization. 

The use of parameterization of point cloud is widespread in computer science 
and engineering. One of the major applications of point cloud parameterization is 
mesh generation. Instead of a convoluted point cloud, mesh reconstruction is usually 
completed on a simpler parameter domain. In mm, Floater and Reimers applied 
their proposed parameterization scheme for meshing point clouds with disk topology. 
In [18], Hormann and Reimers extended the parameterization method in [11] for 
surface reconstruction of point clouds with spherical topology. In [50], Zwicker and 
Gotsmann used their proposed parameterization method for mesh reconstruction of 
genus-0 point clouds. Tewari et al. [44] performed surface reconstruction using their 
proposed doubly-periodic global parameterization. Li et al. m proposed a meshless 
quadrangulation scheme by global parameterization. The input point cloud is cut to 
be with disk topology and parameterized onto the plane for meshing. Zhang et al. [48] 
suggested a mesh reconstruction method of point cloud data by meshless denoising 
and their proposed parameterization scheme. Table |2.1| compares several previous 
works on meshing point clouds using parameterizations. The above previous works 
reflect the importance of parameterization in surface reconstruction of point cloud 
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data. 

Finding a conformal parameterization involves solving differential equations. In 
particular, for conformal parameterizations of point clouds, it is necessary to build 
a discrete analog of the differential operators on point clouds. Numerous works on 
approximating differential operators on point cloud have been reported. In [35], Nay- 
roles et al. described a diffuse approximation method for estimating the derivatives at 
a given set of points. In m, Belkin and Niyogi established a theoretical foundation 
for the Laplace-Beltrami operator on point clouds. Belkin et al. |4] proposed the PCD 
Laplace operator for approximating the LB operator using an integral approximation. 
The moving least square (MLS) method [43l[^ is widely used for the approximation. 

A number of algorithms for the approximation of derivatives are developed based on 
the MLS method [28l |23l |36l [5]. In [63], Lange and Polthier proposed a point set 
analogue of the Laplace-Beltrami and shape operator using the MLS method. In 
[25l l26] , Liang et al approximated the LB operator on point clouds by the MLS 
method with a special weighting function. In [20|, Lai et al. presented a local mesh 
approach for solving PDFs on point clouds. A local mesh structure is constructed 
at each point using local principal component analysis (PCA). Macdonald et al. [33] 
computed react ion-diffusion processes on point clouds. In m, Lozes et al. proposed 
a method to solve PDFs on point clouds for image processing using partial difference 
operators on weighted graphs. 

4. Mathematical background. In this section, we introduce some basic math¬ 
ematical concepts closely related to our work. For more details, readers are referred 

to [4T1 [55] . 

4.1. Conformal maps. An altas of a manifold is said to be eonformal if all of 
its transition maps are biholomorphic. A conformal structure is the maximal con¬ 
formal altas, and a surface with a conformal structure is called a Riemann surfaee. 
Suppose At, M are two Riemann surfaces with local coordinate systems ri(x^, and 
r 2 (x^,x^), where ri,r 2 : ^ are vector-valued functions. The first fundamental 

forms of M and M are respectively defined by 

dslf = Y^gijdx’^dx^, (4.1) 

Consider / : Al -> Af. In local coordi¬ 
nate systems, / can be regarded as / : ^ with f{x^,x‘^) = (/^(x^, /^(x^, x^)). 

The pull-back metric f*ds^ defined on Al, induced by / and is the metric 


ds‘^ 


= 2 ^ 9 ijdx^dx^ and 


M 




where gij = 


_ / dri dri \ ^ _ / dr2 dr2 \ 


dx^ ’ dx^ / ’ \ dx^ ’ dx^ /' 


tdslf = j I dx^dx^- (4-2) 

m,n \ i,j / 

/ is said to be eonformal if there exists a positive scalar function A(x^,x^), called 
the eonformal faetor, such that f'^dsj^ = Ads^. An immediate consequence of the 
above is that every conformal map preserves angles and hence the local geometry of 
the surface. 


4.2. Harmonic maps. By the uniformization theorem, every genus-0 closed 
surface is conformally equivalent to §^. Hence, it is natural to consider mappings 



6 


Choi, Ho and Lui 


between a genus-0 closed surface and the unit sphere. The Dirichlet energy for a map 
f : A4 ^ is defined as 


E{f)= [ (4.3) 

JM 

In the space of mappings, the critical points of E{f) are called harmonie mappings. 
For genus-0 closed surfaces, conformal maps are equivalent to harmonic maps ESI- 
Hence, the problem of finding a conformal map / : ^ is equivalent to an energy 

minimization problem. 

4.3. Quasi-conformal maps. Quasi-conformal maps are a generalization of 
conformal maps. Mathematically, / : C ^ C is a quasi-conformal map if it satisfies 
the Beltrami equation: 


dl 

dz 



(4.4) 


for some complex-valued function ji satisfying ||/i||oo < 1 and ^ is non-vanishing 
almost everywhere. Here, the complex partial derivatives are defined by 


df 1 fdf .df\ , df 1 fdf .df\ 


(4.5) 


/i is called the Beltrami eoeffieient of the quasi-conformal map /. Note that the 
quasi-conformal map / is conformal around a small neighborhood of p if and only if 
li{p) = 0, as Equation (4.4) becomes the Cauchy-Riemann equation in this situation. 

Suppose f : ^ Q 2 and ^ : O 2 ^ ^3 are quasi-conformal maps with the 

Beltrami coefficients /i/ and pig respectively. Then, the Beltrami coefficient of the 
composition map ^ o / : ^ r ^3 is explicitly given by 


M/ + 7;(Mgo/) 
1 + /7/V(Mg ° /) 


(4.6) 


Quasi-conformal maps are also defined between two Riemann surfaces M and N. 
A Beltrami differential on M is an assignment to each chart (I/q,, (fa) of an 

complex-valued function jj^a defined on the local parameter such that iaa{za)j^ = 
lig{zg)^ on the domain also covered by another chart {Ug^'ifg)^ where ^ = -^(fap 
and (fag = 0/3 ° orientation preserving diffeomorphism f : M ^ N is called 

quasi-eonformal associated with ii{z)j^ if for any chart {Ua, (fa) 011 Af and any chart 
{Uglify) on N, the mapping fag := 'ifg o f o f~^ is quasi-conformal associated with 
p(a^^‘ Readers are referred to [14] for more details of quasi-conformal maps. 

4.4. Stereographic projection. In our work, we frequently make use of the 
stereographic projection. Mathematically, the stereographie projeetion is a conformal 
map Pat : ^ C with 


PN{x,y,z) 


X 

l-z 


(4.7) 


The inverse stereographie projeetion is a conformal map Pn" : C ^ with 


PN^{x + iy) = 


2x 


2y 


-1 


1 + + 2/2 ’ 1 + ’ 1 + a;^ + t/^ 


(4.8) 
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Similarly, we define the south-pole stereographic projection Ps : 

X , . y 


Cby 


Ps{x,y,z) = 


1-hz 


-hi 


1 


(4.9) 


The inverse south-pole stereographic projection is the map Pg^ : C ^ with 


PP{x + iy) 


2x 


2 ^ 


1 + ^ 1 + + ^2 


1 — x^ — y^\ 
1 -h x‘^ -h y‘^ J 


(4.10) 


4.5. Point cloud and local system. A point cloud P = 2 ^ 2 , ..., z^} C 

is a set of sample points representing a Riemann surface A4. Because of the absence 
of the connectivity information, we construct a local coordinate system for P on each 
point Zg and approximate the derivatives. To achieve this. We define an atlas (t/g, (js) 
for each point Zg, where Us is an open cover and (js is the associated local coordinate 
function. Us is formed using the collection of all neighboring points of denoted by 
JV{zs). Specifically, we apply the k-Nearest-Neighbors (k-NN) algorithm to define the 
neighborhood. The k-nearest neighborhood M^{zs) of 2 ;^ is a set with the k distinct 
elements in P (including Zs) closest from under the Euclidean 2-norm. In this work, 
we apply the KD-tree implementation by Lin [65] for the computation. We denote 
M^{zs) = { 2 :^, 2 :^,..., z^} with z] = Zg. Then, one common approach for constructing 
a local coordinate system is to define the normal vector as the 2 ;-axis, which is more 
convenient for further computation. There are various methods to obtain the tangent 
planes and the normal vectors for point clouds, such as the principal component 
analysis (PC A) method [T7|. Using the PC A method for 2 ;^, we obtain three vectors 
{e^, e^, e^} which form an ortho normal basis of 

Then, we project J\f^{zs) to the plane spanned by by zl = zl — {zl — 

Zs, el)el, i = 1, 2,..., /c. Now we have the projection M^{zs) = {z\, z^,..., z^} and 
the local coordinates {{xl,yl), {xl^y^),..., {Xg^y^)}, where xl = {zl — Zs^e]) and 
Vs ~ ~ i = 1,2, ...,A:. Therefore, we can define (js : M's ^ by 

(jsi^l) = Also, the neighborhood Af{zs) can be regarded as a graph of its 

projection Af{zs), that is, zl = Zg + xlel + yle^ + fs{xl,yl)el. 


5. An overview of the fast spherical conformal parameterization algo¬ 
rithm for triangular meshes. In this section, we briefly describe the approach in 
[ 6 ] for computing a spherical conformal parameterization of a genus -0 closed trian¬ 
gular mesh M. This approach motivates our proposed parameterization scheme for 

genus -0 point clouds. _ 

To compute a conformal mapping / : M ^ S^, it suffices to solve Equation (4.3). 
This can be achieved by solving the Laplace equation A^/ = 0 subject to ||/|| = 1, 
where A^/ is the tangential component of A/ on the tangent plane of S^. This 
tangential approach was applied by Oberknapp and Polthier in m- Note that this 
problem is nonlinear because of the constraint ||/|| = 1. In [53 Ell, Angenent et al. 
linearize this problem by solving the equation on the complex plane: 


A/ = 0 (5.1) 

given three boundary constraints /(a^) = 6^, where a^, 6^ G C, i = 1,2,3 such that the 
triangle [ai, a 2 , as] and the triangle [5i, 62 , bs] are with the same angle structures. Note 
that A^/ = A/ = 0 since the target domain is now C. As the nonlinear constraint 
II/II = I is removed, the above problem becomes linear and can be solved using the 
cotangent formula m 
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After solving Equation (5.1), the inverse stereographic projection is applied 


for obtaining a spherical parameterization. However, unlike in the continuous case, 
the spherical parameterization in the discrete case is with large conformality distortion 
at the north pole of the sphere due to the discretization and the approximation errors. 
Hence, Choi et al [6] proposed to apply the south-pole stereographic projection Ps to 
map the sphere to a planar domain R C C. Note that the region with large distortion is 
the innermost region of R while the outermost region of R is with negligible distortion. 
Denote the above steps by a map g : M ^ R. To correct the distortion of g, Choi et 
al made use of the quasi-conformal theory. 

Let fig-i be the Beltrami coefficient of the map g~^. Fixing the outermost region 
on i?, Choi et al [6] composed the map g with a quasi-conformal map h : R ^ E>‘^ 
with the associated Beltrami coefficient Let h = u iv and gh = P P 

ir. Specifically, by considering the Beltrami Equation (4.4), each pair of the partial 
derivatives Vx^Vy and Ux^Uy can be expressed as linear combinations of the other m, 


-Vy = aiUx + a2Uy] -Uy = aiVx + a2Vy] 




where ai = ;a 2 




^x — ^2^x T 


l-p‘ 


1 — 


-. Since V 


(5.2) 


= 0 and 


= 0 , the map h can be constructed by solving the following equations 


V- A 


= 0 and V• A 


= 0 


(5.3) 


where A = 


In the discrete case, the above elliptic PDFs (5.3) can 


ai ol2 
0L2 ^3 

be discretized into sparse symmetric positive definite linear systems as described in 
[321 E]. In [58], Jones and Mahadevan derived the system (|5.2|) from the conjugate 


Beltrami equation || = and proposed an alternative approach for solving 

the system. Specifically, the authors considered minimizing the following functional 


^ [ {ai\\/u\‘^ p2a2'Vu'\/v P a3\\/v\‘^) dS, 
^ Jr 


(5.4) 


using a Euler-Lagrange variational approach. Despite the different implementations, 
both of the two abovementioned methods ef fecti vely solve for a quasi-conformal map 
on triangulated meshes. Then, by Equation (4.6), the composition map hog : M ^ 
is with the Beltrami coefficient ghog = 0 and hence h o ^ is conformal. Readers are 
referred to |6] for more details. 

Note that a key step above is the computation of the quasi-conformal map h for 


improving the conformality, which is guaranteed by the composition formula (4.6). 


However, the Beltrami coefficients in the above algorithm are approximated on the 
triangular faces of a mesh. Hence, the above algorithm cannot be directly applied 
for point clouds. Moreover, even if we can define the discrete Beltrami coefficients on 
point clouds. Equation (4.6) may not hold anymore. Therefore, we need to replace 
this key step by a new method suitable for point clouds. 

6 . Meshing genus-0 point clouds using spherical conformal parameteri¬ 
zation. In this section, we discuss our proposed framework for meshing genus-0 point 
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Weight 

Formula of w{d) 

Constant weight 

w{d) = 1 

Exponential weight 

w{d) = exp 

Inverse of squared distance weight 


Wendland weight [46] [47] 


Special weight [25] 

s 

II 

4k II 

0 0 


Table 6.1 

Some common weighting functions for the MLS approximation. 


clouds. The main steps involved include solving a series of Laplace equations on the 
complex plane for the spherical conformal parameterization of a genus-0 point cloud, 
and creating a mesh structure with the aid of the global parameterization. 

6 .1. Approximation of the Laplace-Beltrami operator. In this subsection, 
we explain our approximation scheme for the Laplace-Beltrami operator in the Laplace 


Equation (5.1) on a point cloud P by the moving least-square method. The moving 


least-square method is widely used for approximation [28l [23l |36l O |63l |25l [26] . In 
particular, Liang et al [25l |26] approximated the LB operator on point clouds using 
the MLS method with a special weight function. Our approximation scheme is built 
upon the method in [63] [25] [26] . In this work, we propose a new weight function to 
achieve a more accurate approximation of the LB operator. 

First, we discuss our approximation method for the derivatives on the point cloud 
P = {zi^Z 2 ,... ,Zn}. To simplify the discussion, we only consider the approximation 
on the patch M{zs) of a point 2 ;^ G P. Recall that M{zs) can be regarded as a 
graph of its projection A^(zs), that is, zl = 2 ;^ -\-x\e\ Pyle^ -hDenote the 
derivatives of fs along the -direction and the e^-direction by fsx and fsy respectively. 
We select a set of basic functions , • • •,as a basis and write fs{x,y) ~ 

where {ci}^^ are some coefficients to be determined. In our work, 
we use {1,T, as the basis of the space of all polynomials with second 
order or below, which means m = 6. We add a remark here that m = 6 is an 
appropriate choice for our approximation. Since second derivatives are considered in 
approximating the LB operator, polynomials with at least second order are needed. 
On the other hand, if we fit a polynomial with third order {m = 10) or higher, it will 
be too sensitive to noises and the approximation gets worse. Therefore, m = 6 is a 
suitable dimension for our approximation. 

In the approximation, we aim to minimize 


'^Cjp{xi,yi) - fs{xi,yi) 


i=l 


U = 1 



( 6 . 1 ) 


where Wi = w{\\zi — Zs\\) for some weighting function re : M ^ M. The weight function 
w significantly affects the accuracy and robustness of the approximation. Hence, 


one must carefully choose a suitable weight function. Table 6.1 lists some common 
weighting functions. 
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Note that the information provided by the data points near the center point 2 ;^ 
should be more reliable than that of the data points distant from The closer the 
data points are to Zi^ the more reliable they are. Hence, it is natural to consider a 
smooth weight function which concentrates at . This motivates us to use of a weight 
function of the Gaussian type. As a remark, in Eaii, Belkin et al used a Gaussian 
weight function in the form of exp (—||x — ^|p/4t) for integral approximation. In our 
MLS approximation, we propose another Gaussian-type weight function: 


Ws = 'ie(O) = 1 


Wi = W{ Zi - z. 




Vk . 


for all i ^ 


( 6 . 2 ) 


where h is the maximum distance from Zg in M^{zs). Numerical experiments are 
demonstrated in Section to support our proposed weight function with the specific 
factor ^/k/h‘^ inside the exponent. It can be observed that our proposed weight results 
in more accurate approximations of the LB operator on point clouds. 

With the proposed weight function, we now solve the minimization problem (6.1). 
fl{xi,yi) and fs^ 


Denote /; 


j 

s,i 


Vi)- 



ffh 

fh ■ 

■■ fZA 


fwi 

0 • 

• ■ 0 \ 


/Cl\ 

Let A = 

fl2 

fh ■ 

-pm 
■ ' /s,2 

,D = 

0 

W2 • 

• ■ 0 

, C = 

C2 


[fin 

flu ■ 

. . fm 

J s,n/ 


U 

0 • 



\C772 / 


and 


b = 


ffsA 

fs,2 
\fs,nj 


The minimization problem in (6.1) can be written as follows: 


min ( D(Ac — b),Ac — b 

ceR'^ \ 

We can solve it using the least-square method, namely solving 

D Ac = A^Db. 


(6.3) 


(6.4) 


Next, for any function u defined on the neighborhood M{z)^ we can approximate 
it by a combination oi {fI, 


u = fs{x,y) « '^Cifl{x,y). 


i=l 


Similarly, the coefficients A can be approximated. Let A = 


/wi 

0 • 

■ ■ o\ 


(a\ 


/mi\ 

0 

W2 • 

• • 0 


C2 


U2 

\o 

0 • 

■ • WnJ 

, C = 

\^m / 

, and u = 

\Un) 


D = 


Ci by solving the following least-square problem 

A^DAt = A^Du. 


(6.5) 


Let A = 

(fU 

fh 

fh ■ 
fh ■ 

-pm 
■ ' /s,2 


[fin 

fin ■ 

. . fm j 

J s^n/ 


We can find the coefficients 


(6.6) 
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Since we know the explicit formula of the derivatives of each /], we can compute 
the approximated derivatives of such as 


du ^^dfl , du ^^dfl 


(6.7) 


Now, we are ready to introduce the construction of the LB operator of a smooth 
function u on Af{zs). For any smooth real-valued function u on the Af{z), the LB 
operator of u is given by 


^ W'tl 9i(g"^Wdj{u{z))), ( 6 . 8 ) 

^i=i 

where 2 ) is a point in Af{z)^ {gij) is the metric of the surface at z, VF = ^/det{gij)^ 
and {g 9 ) = ( 5 ^)“^ 

Since zl = fs{x\^y\)) and M{zs) is a graph of JV{zs), we have 


idij) — 




{fs)x{fs)i 


{fs)x{fs)y 

1 + (/.)? 


and 




ifs)l 
ifs)xifs)y 


1 + (L)l -{fs)x{fs)y 

l + (/a)l 

(6.9) 

where W = 

We use Equation ( |6.7| ) to calculate the first order partial derivatives of fg. Then, 
we proceed to compute Au{zs). Since we have a closed form of Au and the LB 
operator is a second order differential operator, by differentiating Equation ( 6 . 8 ), we 
get 


A / X du , . du , . d‘^u , , 

Au{Zs) = ai^{Zs) ^ a2-^[Zs) ^as-^{Zs) 


d‘^u d‘^u 

■ C^4Tr-^_[Zs) + C^s^{Zs) 


dxdy 


' dy‘^ 


( 6 . 10 ) 


where o^i, 0 ^ 2 , trs, 0 ^ 4 , 0^5 are coefficients which depend on partial derivatives of /g. 
This completes our approximation scheme for the LB operator on point clouds. With 
this approximation, we are now ready to describe our proposed spherical conformal 
parameterization algorithm for genus -0 point clouds. 


6 .2. Spherical conformal parameterization of genus-0 point clouds. In 

this subsection, we introduce our proposed method for the spherical conformal pa- 
rameterizations of genus -0 point clouds. 

Given a point cloud P sampled from a genus-0 closed surface At, our goal is to 
find a conformal map / : P ^ which effectively resembles the conformal map 
/ : Ad ^ §^. By the previous section, we can approximate the LB operator A on P. 
Denote the approximated LB operator on the point cloud by Apc- The approximation 
allows us to solve the Laplace equation (5.1) on point clouds for a map 0 : P ^ C. 
More specifically, we solve the following equation 


Apc0 = O (6.11) 

subject to the constraints 0(a^) = hi for i = 1, 2, 3, where ai^hi G C. The choice of the 
three boundary points ai, a 2 , as affects the conformality of the map (j). In the case of 
triangular meshes, ai,a 2 ,a 3 are chosen to be the three vertices of the most regular 
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triangle among all triangles on the input mesh [ 6 ]. Here, the regularity of a triangle 
[ai,a 2 ,a 3 ] is defined by 


Regularity[ai, a 2 , as] = 


TT 

^“3 


+ 


'’-i 


+ 


7 - 


( 6 . 12 ) 


where a, P and 7 are the three angles in the triangle [ai, a 2 , as]. However, in the case 
of point clouds, we do not have the required connectivity information. Hence, we 
choose the three points ai, a 2 , as in a different way. 

Recall that in approximating the LB operator, it is necessary to find the k nearest 


neighboring data points 


, Zg for each point 2 :^ on the point cloud P. We 


consider forming a triple using and two other neighboring points zl and z ^, where 
i j. Different combinations of i and j result in different triples [ 2 ;^, z*, zi]. Then, we 
propose to choose the three boundary points ai,a 2 ,a 3 in the constraint of Equation 
( 6 . 11 ) by considering 


min Regularity [ 2^5 

s,i,j 


7 7^ 


(6.13) 


among all combinations of 5 , i and j. 


After solving Equation (6.11) with our proposed boundary constraints, we apply 


the inverse stereographic projection on 0 (P) to obtain a spherical point cloud. 
Recall that the conformality distortion around the north pole is large due to the 
approximation error in the stereographic projection. Note that the key step in the 
method in [ 6 ] for correcting the distortion via a composition of quasi-conformal maps 
does not work for the case of point clouds. Now, we propose a new method to correct 
the conformality distortion by solely using the LB operator. 

We begin with the south-pole stereographic projection Ps to project the spherical 
point cloud back onto the complex plane. Under the projection, the North pole of 
the sphere, which corresponds to the outermost region of 0(P) C C, is mapped to the 
innermost region on the complex plane. It follows that the outermost region is now 
with very low distortion while the innermost region is with large distortion. We use 
the outermost low-distortion data points as the boundary constraints and solve the 
Laplace equation 7 /; : {Ps o P^^ o (j)){P) C again: 


ApcV^ = 0 


(6.14) 


subject to the boundary constraints 2 p{x) = x for all data points x in the outermost 
low-distortion region. The low-distortion boundary constraints provide us with a more 
accurate result in the inner part of the planar region. Einally, we apply the inverse 
south-pole stereographic projection P^^ and obtain a composition map 


f = O^oPgoPj 


N 


(6.15) 


This step effectively replaces the step in the mesh parameterization algorithm in [ 6 ] 
which involves computing a quasi-conformal map. 


Altogether, by solving Equation (6.11) and Equation (6.14) and using a number 
of projections, we can obtain a conformal map / : P ^ S^. Note that the method 
in [ 6 ] is based on certain manipulations of Beltrami coefficients and quasi-conformal 
maps. In contrast, our new method only involves solving Laplace equations. The 
equivalence between the two approaches can be explained as follows. 

In the first step, the conformality distortion of the spherical parameterization 
is due to the error in the stereographic projection. Then in the approach in [ 6 ], 
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the entire initial parameterization result is used in Equation (5.3) for computing a 


quasi-conformal map in order to cancel the distortion. The method is theoretically 


guaranteed by the composition formula (4.6) of quasi-conformal maps. In contrast, in 


our new approach, we only make use of the most accurate part in the initial parame¬ 
terization result. More explicitly, we use the southern-most regions as the boundary 
constraints and compute the remaining part of the spherical parameterization again, 
with the aid of the LB operator. The replacement of the south-pole step in [6] by our 
new south-pole step can be justified by the following theorem. 

Theorem 6.1. Let {Si,a\dz\‘^) and { 82 ^ p\dw\‘^) be two Riemann surfaees, and p 


is a presented Beltrami differential on Si. Then, the map solved by Equation 
is a harmonie map between {Si, \dz + iidz\^) and {S 2 , p\dw\‘^). Consequently, solving 
the Laplaee equation { 6 A^ is equivalent to solving Equation ( f5.5| ). 

Proof. Let ( be the coordinates of Si with respect to the distorted metric \dz + 
pdz\‘^. Then, the harmonic map between {Si, \dzpdz\‘^) and { 82 ^ p\dw\‘^) is a critical 
point of the following energy 


Eharm{h)= [ p{h{C)){\h^\^ + \h^\'^)dxdy. 
JSi 


(6.16) 


On the other hand, by definition of the Beltrami equation (4.4), the solution to Equa¬ 


tion (5.3) is the critical point of the following energy functional 


EQcif)=[ p{f{z)){\r-Iif^f)dxdy. 

JSi 


(6.17) 


It is shown in m that the above two energy functionals have the same set of critical 
points. Hence, solving the generalized Laplace equation (5.3) for a quasi-conformal 


map is equivalent to solving the Laplace equation (6.14) for a harmonic map under 


the distorted metric. Then, the conformality of our approach is again guaranteed by 


the composition formula (4.6) of quasi-conformal maps. □ 


Therefore, in the continuous case, under suitable boundary conditions in Equation 


(5.3) and Equation (6.14), both of the two methods are theoretically guaranteed for 
producing a conformal map. 

However, in the discrete case, the two methods perform differently. Eor the case 
of triangular meshes, the Beltrami coefficients can be accurately computed and the 


composition formula (4.6) of quasi-conformal maps is accurate under the discretiza¬ 


tion. In this situation, the method in [6] can be effectively applied. Yet, for the case 
of point clouds, we only have an approximation of the LB operator but not the Bel¬ 


trami coefficients, and there is no guarantee about the composition formula (4.6) of 


quasi-conformal maps. Consequently, it is more suitable to use our proposed method 
as it only involves solving the Laplace Equation (6.14). However, since the accuracy 


of our proposed method depends on the boundary constraints in solving the Laplace 
Equation (6.14), the boundary constraints obtained from the initial parameterization 
result may contain small error and hence slightly affect the result in solving Equation 
(6.14). Therefore, it is desirable to perform some more iterations for obtaining a more 


accurate result. 

It is noteworthy that in the parameterization algorithm in [6] for triangular 
meshes, no further steps are required after the second step. However, because of 
the abovementioned issue about the boundary constraints in the Laplace Equation 


(6.14), further iterations are necessary for enhancing the parameterization result. We 


call them the North-South (N-S) reiterations. In each N-S reiteration, two Laplace 
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equations are solved again after the north-pole stereographic projection and the south- 
pole stereographic projection respectively. For solving each Laplace equation, we fix 
the outermost r% points on C to ensure the existence of the solution. 

More specifically, in each N-S reiteration, we first project the previous spherical 
parameterization result onto the complex plane^using the_north-pole stereographic 
projection. Next, we compute a harmonic map 0 : (P/v o /)(P) ^ C by solving the 
Laplace equation 

A,pc(t> = 0 (6.18) 

with the boundary constraints (/)(x) = x for the outermost r% of the data points 
on C. After obtaining 0, the inverse north-pole stereographic projection is again ap¬ 
plied, followed by the south-pole stereogr^phic projection. Then, we compute another 
harmonic map ip : {Pg o py 0(j)0PN o f){P) ^ C by solving the Laplace equation 

= 0 (6.19) 

with the boundary constraints '0(x) = x for the outermost r% of the data points on 
C. We then define the updated spherical parameterization by the composition map 

P“^ o%jj o Ps O P~^ o (j )0 P/v o /. (6.20) 


We check whether the a^ove updated parameterization result is close to the previous 
parameterization result /. If yes, then the parameterization is stable and we complete 
the algorithm. If no, we apply another N-S reiteration on the updated parameteriza¬ 
tion point by repeating the procedures and so on. In practice, we choose r = 10. Our 
proposed spherical conformal parameterization scheme is outlined in Algorithmic 
To explain the motivation of our proposed N-S reiteration scheme, we define the 
N-S Dirichlet energy by 

E{f) = i (E(Piv(/)) + EiPsim , (6.21) 


where P(/) is the Dirichlet energy. It follows that the minimum of E is attained if and 
only if P(P/v(/)) and E{Ps{f)) are minimized, which implies that E{f) is minimized 
and / is conformal. Therefore, to find a conformal /, we can consider minimizing the 
N-S Dirichlet energy E{f). More specifically, we aim to minimize both P(P/v(/)) and 
E{Ps{f)). Note that these two energies are respectively minimized by solving the 
Laplace equations ( |6.18 ) and (6.19) in our proposed N-S reiteration. Introducing the 
N-S reiteration for minimizing the energies is advantageous for two reasons. Firstly, 
it linearizes the computation as we only need to solve linear systems on C. Secondly, 
it avoids the error induced by the stereographic projection as we consider both the 
north-pole step and the south-pole step in each reiteration. 

Note that theoretically we only need to take away the two points of infinity, and at 
least 3 points on C are required to be fixed to guarantee the existence of the solution 
of the two Laplace equations. However, in terms of the numerical computations, 
the large matrix equations may be ill-posed if only 3 points are fixed as boundary 
constraints. Therefore, we consider fixing the outermost r% of data points in solving 
Equations (6.18) and (6.19). Nevertheless, fixing these extra points may not affect 
the accuracy of the solution too much. This can be explained with the aid of the 
Beltrami holomorphic flow [14] as follows. 

Theorem 6.2 (Beltrami holomorphic flow on C [E]). There is a one-to-one eor- 
respondenee between the set of quasi-eonformal diffeomorphisms ofC that fix the points 
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0,1 and oo and the set of smooth eomplex-valued funetions on C with 


k <1. 


Furthermore, the solution to the Beltrami equation depends holomorphieally 

on fi. Let be a family of Beltrami eoeffieients depending on a real or eomplex 

parameter t. Suppose also that ii{t) ean he written in the form 


li{t){z) = ii{z) + tu{z) + te{t){z) (6.22) 

for z G C, with suitable /i in the unit ball o/(7^(C), z/, e(t) G I/^(C) sueh that 
||e(t)||oo ^0 as t ^ 0. Then, for all w G C, 


/A‘(*)(^) = + tV{r, u){w) + 0 {\t\) 

loeally uniformly on C as t ^ 0, where 

-1) f v{z){{f^^)^{z)r 


(6.23) 


Vir,u)iw) = --- 


L 




dxdy. 


(6.24) 


In our case, as the conformality distortion of the outermost region is negligible, u 
is compactly supported around origin. Hence, it can be deduced from Equation (6.24) 
in the above theorem that the data points located farther away from the origin will be 
associated with a smaller flow V, as the denominator in the integral becomes larger. 
Therefore, in each iteration, the outermost points will remain almost unchanged, while 
the innermost points (which have the largest conformality distortion) will be adjusted 
and improved. In other words, fixing more outermost points for the numerical stability 
does not affect the solution much. Numerical experiments are presented in Section 
[^to verify the convergence of our N-S reiteration scheme. Intuitively, the boundary 
constraints in Equations (6.18) and (6.19) are adjusted to the positions associated 
with a conformal map by the iterations. They are observed to eventually stabilize 
and hence we obtain the desired conformal map by solving Equations (6.18) and 
(6.19) with these boundary constraints. 

Einally, we make an important remark about our proposed spherical conformal pa¬ 
rameterization algorithm for genus-0 point clouds. In addition to genus-0 point clouds, 
our proposed algorithm also efficiently works on genus-0 triangular meshes. Note that 
for triangular meshes, the LB operator can be easily constructed by computing the 
cotangent weights on the given mesh structures. Also, solving Laplace equations on 
the complex plane requires only linear time. Hence, our proposed algorithm can serve 
as an alternative approach for computing spherical conformal parameterizations of 
genus-0 closed triangular meshes in linear time. 


6.3. Improving the distribution of the spherical parameterization. It 

is obvious that spherical conformal parameterizations are unique only up to Mobius 
transformations. Although the conformality does not change under the Mobius trans¬ 
formations, the distribution of the points on the sphere does. The distribution is 
crucial for meshing. Hence, it is desirable to obtain an even distribution of the points 
on the sphere. 

In the spherical conformal parameterization algorithm for triangular meshes [6], 
Choi et al proved the following theorem: 

Theorem 6.3 (See [6], P.75). Let Ti and T 2 be two triangles ofC. The produet 
of the perimeters of Ti and Ps{Pf^^{T 2 )) is invariant under arbitrary sealing of Ti 
and T 2 . 
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Algorithm 1: Our proposed spherical conformal parameterization algorithm. 
Input: A genus-0 point cloud P. 

Output: A spherical conformal parameterization / : P — 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


Approximate the LB operator on P and denote it by Apc; 


Find the most regular triple of points by solving problem (6.13); 
Obtain a map 0 : P ^ C by solving the Laplace equation (6.11); 
Apply the inverse stereographic projection P^^ : C 
Apply the south-pole stereographic projection Ps : 


on 0(P); 
^ C on {P^^ 


^){py. 


Solve the Laplace equation (6.14) for -0 : {Ps o P~ o (j)){P) C; 

Apply the inverse south-pole stereographic projection P^^ and denote the 
overall composition of the maps by / = Pg^ o'll; o Ps oPn^ 

repeat 

Update / by /; 


o0; 


Solve the Laplace equation (6.18) for (p : (P/v o /)(P) ^ C; 

Solve the Laplace equation (6.19) for ip : {Ps o P^^ o p o Pjs^ o /)(P) ^ C; 


Update / by P^ o -0 o P^ o P^^ o p o P/v o /; 
uutil mean{\\f{pi) - f{pi)\\‘^) < e; 


With this theorem, Choi et al achieved an even distribution of a spherical param¬ 
eterization mesh by applying the stereographic projection on the sphere, and then 
considering the outermost triangle T and the innermost triangle t on the complex 
plane. They scaled the planar domain by a factor so that T and t are with the same 
perimeters on the sphere, under the inverse stereographic projection. 

In our case, the above idea does not work as we do not have any information 


about the connectivity of the point clouds. However, we can extend Theorem 6.3 for 
point clouds by considering two sets of points. The extension is as follows: 

Theorem 6.4. Let {ui]^^ and {vj}'j^Q be two sets of points on C. Then 


Y, l|A«i - Atioll E ||Ps(P^'(At;,)) - P5 (P^'(A^;o))|| 


U = 1 


= ( \\Ui - noil ) I \\Ps{Ppivj)) - PsiPi^\vo))\\ 

Vi=i 




for any sealing faetor A 7^ 0. In other words, the produet is an invarianee under 
arbitrary sealing. 


Proof We prove the theorem using the approach in [6]. Note that for any 2: = 
X + iy, we have 


PsiPpiz)) = Ps{Pi^\x + iy)) 

_ 2x 2y 

_ 1 +x^+y^ . l+x‘^+y‘^ 

-j I -l+x‘^+y‘^ 1 I -l+x‘^+y‘^ 

-L l+x^+y-2 -L l+x^+y-2 

^ -X . y ^ -Re{z) . Im{z) 

x2 ^ y2 * 3 , 2 _|_y 2 


(6.25) 
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Hence, for any scaling factor A 7 ^ 0, we have 


^ ||A«, - Attoll I I ^ - PsiPi^\Xvo))\\ 

i=i 


j=i 


= (y]iiAMi-AMoii) 


a=i 


u=i 


—Re{Xvj) ^ .Im{Xvj) —Re{Xvo) ^ .Im{Xvo) 


lAtijP lA^iP |Atiop |Atio|" 


= (a^ IK -Moll 1 I 


i=l 


i=l 


—Re{vj) ^ .Im{vj) —Re{vo) ^ .Im{vo) 


i.|2 


i .|2 


l^ol^ 


It'ol^ 


= ( ^ ||m, - noil ) I ||Ps(P^iK)) - Ps(P^'K))|| I . 

U = 1 


(6.26) 


Therefore, the product is an invariance. □ 

To apply this theorem for obtaining an even distribution of our spherical parameter¬ 
ization result, we propose to use the average distance between the poles on the unit 
sphere and their k-NN neighborhoods. More specifically, suppose vn and vs are the 
northernmost point and the southernmost point on the spherical parameterization 
result f{P) obtained by Algorithm respectively. By the north-pole stereographic 
projection P/v, is mapped to the point xn on the complex plane. On the other 
hand, by the south-pole stereographic projection P/v, vs is mapped to the point xs 
on the complex plane. Denote the average distances of xn and xs to their k-NN 
neighborhood on their corresponding planar domain by and ds respectively, djsf 
and ds are explicitly given by 


dp = mean({|Pjv(/(^)) -xn\-. z G ^(vjv))}) 

and 

ds = mean({|Ps(/( 2 ;)) -xs\-. z G 

Then, we scale the whole planar domain (P/v o /)(P) by a scaling factor 

^dp X dg 


(6.27) 


(6.28) 


(6.29) 


Now, denote the two updated average distances by dp and dg. It follows that 


\/dp X ds , - 

dp = Xdp = — - - X dp = ^ dp X dg. 

dp 

(6.30) 

Also, by Theorem |6.4[ we have 


dp X ds = dp X dg. 

(6.31) 

Therefore, 



dg = dp X dg X — = 'JdpX dg. 
dp 


(6.32) 
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In other words, the two updated average distances dp and ds defined on the new 
spherical parameterization result Pn^ (A(P/v(/(P)))) are equal. This indicates that 
the distribution of the points at the two poles is balanced. Hence, Algorithm 
together with the described balancing scheme provide us with a a spherical conformal 
parameterization with an even distribution. Our balancing scheme is summarized in 
Algorithm 


Algorithm 2: Our proposed balancing scheme for better distribution. 

Input: A spherical conformal parameterization / : P ^ §^. 

Output: A spherical conformal parameterization with improved distribution. 

1 Apply the north-pole stereographic projection Pn on /(P); 

2 Denote the northernmost and the southernmost points of /(P) hy vn and vs- 

Multiply all points in Pjv(/(P)) by a scaling factor A = where 

dp = niean({|Pjv(/(2)) - xn\ ■ z e V'‘(/“^(vjv))}) and 
ds = mean({|Ps(/(z)) - xs\ : z e J^^if~^ivs))}y, 

3 Apply the inverse north-pole stereographic projection Py on A(Piv(/(P))); 


6.4. Meshing using spherical conformal parameterization. In this sub¬ 
section, we present our meshing framework for genus-0 point clouds. Directly trian¬ 
gulating a point cloud is difficult because of its complicated geometry. However, with 
the aid of the spherical conformal parameterization of point clouds, the difficulty is 
significantly alleviated. Instead of triangulating a point cloud, we triangulate the unit 
sphere obtained by our spherical conformal parameterization algorithm. Algorithms 
for triangulating a spherical point cloud are well-established. In particular, the spher¬ 
ical Delaunay triangulation method, which computes a Delaunay triangulation on the 
unit sphere, is the most suitable method for our purpose. 

Delaunay triangulations are widely used in computer graphics because of their 
good triangle quality. More specifically, Delaunay triangulations are advantageous as 
they maximize the minimum angle in every triangle in the triangulations and hence 
avoid skinny triangles. With this important property, the triangulations generated by 
this method are more regular than the common triangulation methods. 

By applying the spherical Delaunay triangulation method on the spherical con¬ 
formal parameterization of a genus-0 point cloud, we obtain a nice triangulation on 
the spherical point cloud. Since the points on the original point cloud and those 
obtained by the spherical conformal parameterization have a I-1 correspondence, the 
triangulation on the spherical point cloud naturally induces a triangulation on the 
original point cloud. It is noteworthy that since the parameterization is conformal, 
the angles of the new triangulation on the original point cloud are well preserved. In 
other words, the regularity of the triangulation defined on the original point cloud 
closely resembles that of the spherical Delaunay triangulation. Moreover, the meshing 
result is guaranteed to be a genus-0 closed triangular mesh because of the spherical 
Delaunay method. This completes our goal of meshing a genus-0 point cloud. Our 
meshing framework is described in Algorithmic 

Before ending this section, we make an important remark about a possible exten¬ 
sion of our proposed framework. In fact, our proposed parameterization and meshing 
scheme can be possibly extended for point clouds with disk topology. In this case, 
we can first extend the double covering technique [HI 13 to turn a point cloud with 
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Algorithm 3: Our proposed meshing framework for genus-0 point clouds. 
Input: A genus-0 point cloud P. 

Output: A triangular mesh M = (P^T) where T is a triangulation of P. 

1 Apply Algorithm and Algorithm to obtain a spherical conformal 
parameterization / : P ^ 

2 Compute a triangulation T on /(P) using the spherical Delaunay algorithm; 

3 Use T to form a triangular mesh M = (P^T); 


disk topology into a genus-0 point cloud. More specifically, given a point cloud P of 
a simply-connected open surface Al, we turn P into a point cloud P with spherical 
topology and approximate the derivatives on it by the following steps. 

Step 1: Approximate the derivatives on P using the /c-NN algorithm and the moving 
least square method. 

Step 2: Duplicate P and denote the copy of it by P'. 

Step 3: Define the derivatives on P' using the results in Step 1, with reversed orien¬ 
tations. _ 

Step 4: Identify the boundary point^of P, P' and obtain a genus-0 point cloud P. 
Step 5: Create the LB operator for P using the derivatives on P and P'. 

Then, we^an apply the abovementioned spherical conformal parameterization 
algorithm on P to obtain a spherical point cloud. After that, by applying the stere¬ 
ographic projection on the southern hemisphere, we obtain a planar conformal pa¬ 
rameterization of P. Finally, we can easily compute a Delaunay triangulation on the 
planar parameter domain. Since both the parameterization algorithm and the stereo¬ 
graphic projection produce conformal results, this triangulation on the planar domain 
accurately induces a regular mesh structure on P. This completes the task of meshing 
a point cloud with disk topology. 

7. Experimental results. In this section, we demonstrate the effectiveness of 
our proposed framework for meshing genus-0 point clouds using spherical conformal 
parameterization. In the following, we assess the performance of our proposed frame¬ 
work in different aspects. The datasets used in the experiments are freely adapted 
from the AIM@SHAPE Shape Repository [51], the Stanford 3D Scanning Repository 
[52] and the RGB-D Scenes Dataset v.2 [59|. The mentioned algorithms are imple¬ 
mented in MATLAB. The sparse linear systems for the Laplace equations are solved 
using the built-in backslash operator (\) in MATLAB. All experiments are performed 
on a PC with an Intel(R) Core(TM) i5-3470 CPU @3.20 CHz processor and 8.00 CB 
RAM. 

7.1. The performance of our approximation of the Laplace-Beltrami 
operator. In this work, we apply the MLS method with a new weight function 
for approximating the Laplace-Beltrami operator. It is natural to ask whether our 
proposed weight function produces better results. It is also necessary to compare other 
approximation approaches such as the local mesh method to justify our choice. In 
this subsection, we compare the numerical accuracy of the local mesh method and the 
MLS method with several weighting functions for approximating the LB operator on 
point clouds. More specifically, we compare the performance of the following methods: 

1. The local mesh method [20] . 

2. The MLS method with the Wendland weight in @61117], 

3. The MLS method with the Gaussian weight in Ezna, 






Fig. 7 . 1 . Two experiments for assessing the approximation aecuracies of the Laplaee-Beltrami 
operator. In eaeh experiment, we generate a point eloud on the unit disk and transform it using a 
conformal map with an explicit formula. We then approximate the LB operator on the transformed 
point cloud and solve the Laplace equation hack onto the unit disk. Top: the first experiment. 
Bottom: the second experiment. 


4. The MLS method with the special weight in m, and 

5. The MLS method with our proposed weight function. 

Experiments are carried out for assessing the numerical accuracies of the above- 
mentioned approaches. Figure |7.1| shows the setups in two of the experiments. In 
each experiment, we first generate a point cloud on the unit disk. This serves as the 
ground truth in our analysis. Then, we transform the point cloud using a conformal 
map with an explicit formula. We apply the mentioned approximation schemes for 
approximating the LB operator on the transformed point cloud. Then, we solve the 
Laplace equation with the circular boundary constraints on the original unit disk. 
Theoretically, the result obtained by the disk harmonic map should be exactly the 
same as the original point cloud, as the transformation is given by a conformal map 
with an explicit formula. In other words, the ideal position error between the disk 
harmonic map and the original point cloud should be 0. By measuring the maximum 
and average position error between the pairs of points, we can evaluate the accuracy 
of the aforementioned approximation schemes for approximating the LB operator. 

Table |7.1| illustrates the approximation error of different approaches in the two 
experiments. It is noteworthy that in both experiments, the MLS method with our 
proposed weight function produces approximations which are much more accurate 
those produced by the local mesh method and the MLS method with the Wendland 
weight [46l [47] and the Gaussian weight in isaii]. With similar and negligible average 
position errors, our proposed scheme reduces the maximum position errors by about 
25% on average when compared with the MLS method with the special weight [25] . 
The comparisons reflect the advantage of our proposed method for approximating the 
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Method 

maximum 
position error 

average position 

error 

Local mesh method [20] 

1.3427 

0.0179 

MLS with the Wendland weight in [33 [47] 

3.3074 

0.1696 

MLS with the Gaussian weight in [571 H] 

0.5697 

0.0114 

MLS with the special weight in [25] 

0.0427 

0.0006 

MLS with our proposed weight 

0.0245 

0.0004 


Method 

maximum 
position error 

average position 

error 

Local mesh method [20] 

1.5148 

0.0271 

MLS with the Wendland weight in [46] [47] 

2.0082 

0.0803 

MLS with the Gaussian weight in [571 S] 

1.5460 

0.0925 

MLS with the special weight in [25] 

0.0110 

0.0001 

MLS with our proposed weight 

0.0103 

0.0002 


Table 7.1 

The approximation error in the two experiments. Top: the first experiment. Bottom: the 
second experiment. 


LB operator. 


7.2. Performance of our proposed spherical conformal parameteriza¬ 
tion. After demonstrating the advantage of our approximation scheme for the LB 
operator, we investigate the performance of our proposed spherical conformal param¬ 
eterization algorithm for genus-0 point clouds. Figure |7.2| and Figure 7.3 show the 


results of parameterizing a lion point cloud and a bulldog point cloud using our pro¬ 
posed parameterization method. Two more convoluted examples are shown in Figure 
|7.4| The experiments demonstrate the effectiveness of our proposed algorithm for 
convoluted point cloud data. 

Moreover, with the aid of the spherical conformal parameterization, we can create 
a Delaunay triangulation on the spherical parameterization result by the spherical 
Delaunay algorithm and define an induced triangulation on the input point cloud. 
Using the mesh structures, we can measure the angle differences of the two meshed 
point clouds and hence effectively evaluate the conformality of our parameterization 
scheme. In particular, we define the conformality distortion of the parameterization 
by the angular distortion between the two meshes. The angle difference provides 
an accurate measurement of the conformality distortion of the parameterizations. It 
can be easily observed in Figure [7^ and Figure [73] that the histograms of the angle 
differences highly concentrate at 0. Besides, for better visualizations of the spherical 
conformal parameterization results, we color the surfaces by the approximated mean 
curvature on the source surfaces. It can be observed from the colored figures that the 
local geometries of the point clouds are well preserved under our proposed spherical 
conformal parameterization algorithm. Figureshows the difference plots of several 
experiments. The plots demonstrate the convergence of our North-South reiteration 
scheme. 

We compare our proposed spherical conformal parameterization method for genus- 
0 point clouds with the spherical embedding method proposed by Zwicker and Gots- 
man m and the global conformal map m- In our experiment, k = 2b nearest 
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Point clouds 

No. of 
points 

Performance 

Our 

proposed 

method 

Spherical 

embedding 

m 

Global 
conformal 
map [25] 



Time (s) 

8.1768 

23.7988 

16.7058 

Soda Can 

6838 

Mean( (5) 

0.5902 

4.0431 

2.3352 



SD(|<5|) 

0.8007 

5.2731 

1.9803 



Time (s) 

13.0919 

37.4124 

18.7151 

Hippocampus 

10242 

Mean( ) 

1.2855 

14.3072 

1.3062 



SD(|<5|) 

1.4701 

19.6461 

1.5100 



Time (s) 

30.7785 

87.0887 

40.4214 

Max Planck 

21530 

Mean( ) 

0.7326 

8.6058 

1.0792 



SD(|5|) 

1.0803 

14.0857 

1.5756 



Time (s) 

50.7390 

132.1765 


Cereal Box 

33061 

Mean( (5) 

0.6523 

12.3573 

Fail 



SD(|<5|) 

0.9165 

14.0440 




Time (s) 

114.7057 

291.8312 

122.8818 

Spiral 

48271 

Mean( ) 

0.8580 

16.4704 

0.9658 



SD(|<5|) 

1.3280 

22.5073 

1.3135 



Time (s) 

115.5802 

198.2285 

126.7621 

Brain 

48487 

Mean( (5) 

1.4266 

35.2629 

2.2495 



SD(|5|) 

2.9093 

35.7986 

2.7430 



Time (s) 

88.9297 

206.9920 

113.4447 

Bulldog 

49797 

Mean( (5) 

1.5432 

16.2010 

1.8700 



SD(|<5|) 

2.9183 

21.1544 

3.1891 



Time (s) 

95.8935 

212.5685 

136.3296 

Chinese Lion 

50002 

Mean( ) 

1.8474 

19.1579 

2.4907 



SD(|<5|) 

1.9286 

22.7259 

2.6207 



Time (s) 

198.6064 

360.7178 

227.0290 

Bimba 

74764 

Mean( (5) 

0.6227 

18.0340 

0.6379 



SD(|5|) 

0.8129 

20.6272 

0.7975 



Time (s) 

427.7658 

731.8661 

560.6077 

Igea 

134345 

Mean( (5) 

0.7076 

5.0853 

3.8293 



SD(|<5|) 

1.4273 

8.2623 

2.9703 



Time (s) 

676.4106 

995.7537 


Armadillo 

172974 

Mean( ) 

1.4167 

23.2354 

Fail 



SD(|<5|) 

1.6855 

23.9892 




Time (s) 

1305.9013 

1484.7682 

1642.9208 

Lion Vase 

256094 

Mean( (5) 

2.0920 

17.8501 

3.6696 



SD(|5|) 

4.1052 

21.9588 

5.8502 


Table 7.2 

Performances of three spherical conformal parameterization methods for genus-0 point clouds. 
To quantitatively evaluate the conformality of the parameterization, we build a mesh structure on the 
spherical parameterization using the spherical Delaunay method and then create an induced mesh 
structure on the original point cloud. The conformality distortion (denoted by S) of the parameter¬ 
ization is assessed using the angle difference (in degrees) between an angle on a meshed point cloud 
and the mapped angle on the meshed spherical parameterization result. 
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Conformality distortion 



Fig. 7.2. Parameterizing a lion point eloud. Top left: A lion point eloud. Top middle: The 
spherical conformal parameterization of the lion point cloud. Top right:Middle left: A triangulation 
created by our method. Bottom left and middle: The triangulated point cloud and the spherical 
parameterization result colored with the approximated mean curvature at each vertex. Bottom right: 
The conformality distortion of the parameterization based on the triangulation. 


neighbors of every point are used for approximating the LB operator in Algorithm 
The stopping threshold for the N-S reiteration is set to be e = 0.0001. Table 
summarizes the computational time and the conformality distortion of the three 
schemes. It is noteworthy that our proposed method produces spherical conformal 
parameterizations with better conformality. The better conformality obtained by our 
method is attributed to the south-pole step in our algorithm, which is conceptually 
equivalent to a composition of quasi-conformal maps. With this specific step, our 
method can hence further reduce the conformality distortion of the parameterization. 
Moreover, our method is more efficient than the algorithms in m and [25] . The above 
results indicate that our parameterization algorithm preserves the local geometry of 
the point clouds very well. 



7.3. Performance of our meshing scheme. As mentioned in the last sub¬ 
section, we generate mesh structures on genus-0 point clouds by building Delaunay 
triangulations on their spherical conformal parameterizations. Our meshing scheme 
has two important advantages. First, the regularity of the triangulations generated 
is guaranteed by the preservation of the angle structures of the Delaunay triangu¬ 
lations computed on the spherical parameterizations. As the angle structures are 
well retained under the spherical conformal parameterization, a regular triangulation 
defined on the parameterized point clouds can effectively induce a regular and almost- 
Delaunay triangulation on the original point clouds. Besides Figure [7^ Figure [73] 
and Figure |7.4[ some more examples of triangulations created by our approach are 
shown in Figure 7.6 It can be observed that our meshing method can handle point 
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Fig. 7.3. Parameterizing a bulldog point eloud. Top left: A bulldog point eloud. Top middle: 
The spherical conformal parameterization of the bulldog point cloud. Top right: A triangulation 
created by our method. Bottom left and middle: The triangulated point cloud and the spherical 
parameterization result colored with the approximated mean curvature at each vertex. Bottom right: 
The conformality distortion of the parameterization based on the triangulation. 


clouds with different geometry. High quality triangulations can be created even with 
the presence of sharp, non-convex and convoluted regions on the input point clouds. 
Second, unlike most of the existing meshing methods, the meshes produced by our 
proposed scheme are guaranteed to be genus-0 closed meshes. No holes or unwanted 
boundaries will be present in our meshing result. Hence, post-processing steps are 
not required in our meshing scheme. 

We compare our meshing method with three existing meshing approaches. As an 
example of parameterization-based approaches, Zwicker and Gotsman m generate 
triangulations for a genus-0 point cloud with the aid of the spherical embedding 
algorithm and the spherical Delaunay triangulation method. On the other hand, 
two typical methods for meshing without using parameterizations are the marching 
cubes algorithm [30] and the Tight Cocone algorithm |8]. Figure |7.7| provides a 
comparison between our method and the three mentioned approaches. It can be 
observed that our meshing scheme and the Tight Cocone algorithm |8] produce high 
quality triangulations, while the triangulations produced by the approaches in m 
and [30] consist of certain sharp and irregular triangles. Also, the result by the 
marching cubes algorithm contains holes while our method is topology preserving. 
Therefore, unlike the marching cubes algorithm, no further post-processing is needed 
in our meshing scheme. 

To quantitatively assess the “almost-Delaunay” property of our meshing results 
induced by the spherical Delaunay triangulations, we recall that a Delaunay trian¬ 
gulation satisfies the opposite angle sum property: For every edge in a Delaunay 
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Fig. 7.4. Parameterizing two eonvoluted point elouds of a human brain and a spiral. Left: The 
input point elouds. Middle: The spherieal eonformal parameterizations obtained by our proposed 
algorithm. Right: The triangulations created by our method. 





Fig. 7.5. Several plots of the difference mean{\\f(pi) — f{pi)\\‘^) with the number of iterations 
in our North-South reiteration scheme. Left: Cereal box. Middle: Hippocampus. Right: Bulldog. 
The plots demonstrate the convergence of our proposed algorithm. 


triangulation, the two angles a and P opposite to the edge satisfy 

(a + /3 < TT. (7.1) 

For the abovementioned meshing approaches, we define the Delaunay ratio of the 
resulting triangulation by 

# of edges with opposite angles a, P s.t. a + /d < tt {12) 

Total of edges in the triangulation 

A higher Delaunay ratio indicates that the triangulation is closer to a perfect Delaunay 
triangulation. The Delaunay ratios by different meshing algorithms are presented in 
Table [73| Because of the high accuracy of our spherical conformal parameterizations, 
our proposed scheme achieves significantly better triangulation results when compared 
with the spherical embedding algorithm [50] and the marching cubes algorithm [30] . 
Also, our results are comparable to or even slightly better than those of the Tight 
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Fig. 7.6. Meshes generated by our proposed method and a zoom-in of them. The regularity 
of the triangulations is attributed to our spherieal eonformal parameterization and the spherieal 
Delaunay method. 


Cocone algorithm [8]. The comparisons demonstrate the advantages of our proposed 
meshing scheme. A further comparison between our method and the Tight Cocone 
algorithm [8] is given in the following subsection. 

In addition, we can generate quadrangulations of point clouds with the aid of the 
spherical conformal parameterization. Two examples of the quad meshes generated by 
our method are given in Figure [^8] To create quad meshes of point clouds, we make 
use of a standard spherical quad mesh and our spherical conformal parameterization 
results. With the aid of the spherical conformal parameterizations, we can interpolate 
the standard quad mesh onto the input point clouds and thus generate quad mesh 
representations. Because of the conformality of our parameterization scheme, the 
resulting quad meshes are with high quality. Also, the meshes are guaranteed to be 
topology preserving. 

Before ending this subsection, we demonstrate the significance of our proposed 
balancing scheme in the spherical conformal parameterization. The redistribution is 
vital for the meshing quality. Figure [73| shows the meshing results with and without 
the redistribution scheme. It can be easily observed that if the spherical parame¬ 
terization of a genus-0 point cloud is unbalanced, then on the mesh generated by 
interpolation with the aid of the spherical conformal parameterization, most of the 
vertices will be concentrated at one small region of the mesh. As a result, most 
features of the underlying surface are lost. In contrast, with our proposed balancing 
scheme, a high quality mesh can be effectively generated. The above results reflect the 
importance of our balancing scheme in the point cloud parameterizations for meshing. 
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Point clouds 

Our 

proposed 

method 

Spherical 

embedding 

m 

Marching 
cubes [30] 

Tight 
Cocone [8| 

Soda Can 

0.99 

0.93 

0.82 

0.98 

Hippocampus 

0.99 

0.91 

0.82 

0.99 

Max Planck 

0.99 

0.93 

0.82 

0.99 

Cereal Box 

0.99 

0.88 

0.81 

0.99 

Spiral 

1.00 

0.85 

0.82 

0.99 

Brain 

0.99 

0.82 

0.82 

0.99 

Bulldog 

1.00 

0.86 

0.82 

0.99 

Chinese Lion 

0.99 

0.84 

0.82 

0.99 

Bimba 

1.00 

0.88 

0.81 

1.00 

Igea 

0.97 

0.90 

0.83 

0.97 

Armadillo 

0.98 

0.80 

0.82 

0.98 

Lion Vase 

0.99 

0.85 

0.83 

0.99 


Table 7.3 

The Delaunay ratios of different meshing approaehes. The ratio assesses the proportion of 
edges in the resulting triangulations that satisfy the opposite angle sum property a ^ < tt in a 
triangulation. A Delaunay ratio exactly equals 1 indicates that the triangulation is Delaunay. 


7.4. Stability under geometrical and topological noises. Our meshing 
framework is stable under geometrical and topological noises of the input genus-0 
point clouds. In some situations, the point clouds obtained by 3D cameras are geo¬ 
metrically noisy. To compute triangulations which represent the underlying surfaces, 
we can first apply a Poisson filtering on the noisy point clouds. Then, with the aid 
of our spherical conformal parameterization, we can obtain high quality triangula¬ 
tions on a uniform spherical point cloud and interpolate them back onto the filtered 
point clouds to produce meshed surfaces. We first demonstrate the effectiveness of 
our algorithm by two synthetic examples. 

Figure 7.10 shows a synthetic point cloud with 3% uniformly distributed random 
noise and our meshing result. We can also construct a faithful triangulated mesh on 
the geometrically noisy point cloud without any filtering or sampling procedure. Fig¬ 
ure |7.11| shows the triangulation result of our meshing scheme and the Tight Cocone 
algorithm [8] on the noisy point cloud in Figure 7.10 All points of the point cloud 


are considered and fixed in the construction of the triangulation. It can be observed 
that there are irregular triangulations and topological holes on the result by [8] , while 
our meshing scheme guarantees a regular and topology preserving triangulation even 
for noisy input point clouds. 

Besides, it is common that the sampling processes result in non-uniformly sampled 
point clouds. In particular, there may be large holes on certain parts of the point 
clouds sampled from genus-0 objects, which create topological ambiguities and hinder 
mesh generations. Our parameterization and meshing scheme produce satisfactory 
results with these topological noises. Moreover, the meshes generated are guaranteed 
to be genus-0 closed meshes. Figure 7.12 shows a synthetic point cloud with 1021 
randomly created topological holes. It can be observed that our algorithm produces 
a satisfactory meshing result. 

Then, we apply our algorithm for real 3D scanned noisy point cloud data. Several 
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Fig. 7.7. A comparison between our meshing seheme and other approaehes. A front view of 
the triangulated point eloud and a zoom-in of the nose are shown for each method. Left to right: 
Our meshing result, the method in m, the marehing eubes algorithm and the Tight Coeone 
algorithm 



Fig. 7.8. Quad mesh generation on point elouds using our proposed method. 


raw point cloud data are adapted from the RGB-D Scenes Dataset v.2 [59]. Figure 


The above experiments demonstrate the stability and robustness of our proposed 
method for noisy point clouds. 

7.5. Multilevel representations of genus-0 point clouds. With our pro¬ 
posed spherical conformal parameterization scheme, multilevel representations of a 
genus-0 point cloud can be easily achieved. We start with a coarse spherical point 
cloud. The vertices on the sphere can be interpolated onto the genus-0 point cloud 


7.13 shows two point clouds of a soda can and a cereal box, and our meshing results. 
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Fig. 7 . 9 . The effect of our balancing scheme on meshing a genus-0 point cloud of Max Planck. 
Top left: A spherical conformal parameterization without the balancing scheme. Bottom left: A 
spherical conformal parameterization with the balancing scheme. Middle: The front view of the the 
meshing results by interpolation with the aid of the parameterizations. Right: The back view. 


with the aid of its spherical parameterization. Then, we can progressively subdi¬ 
vide the sphere using existing subdivision methods, such as the butterfly subdivision 
method [9] and the loop subdivision method [29] . For each subdivided sphere, we can 
repeat the mentioned interpolation procedure and obtain a coarse representation of 
the point cloud. This method results in multilevel representations of the point cloud. 
As the subdivision level increases, more details of the point cloud are represented. 
Examples of multilevel representations of genus-0 point clouds are given in Figure 


subdivision method [29|. The subdivision connectivity of the results can be easily 
observed. The results indicate that our method can effectively generate the multilevel 
representations of genus-0 point clouds. 

8. Conclusion and Future Work. In this paper, we presented a novel frame¬ 
work for meshing genus-0 point clouds via global spherical conformal parameteri¬ 
zations. We extended and improved the parameterization algorithm for triangular 
meshes in [6]. Firstly, we enhanced the accuracy for approximating the LB operator 
on point clouds by using a new Gaussian-type weight function. Secondly, we pro¬ 
posed an iterative scheme called the N-S reiteration to replace the step of solving 


7.14|and Figure [7T^ In our examples, the subdivisions are generated using the loop 
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Fig. 7.10. Meshing a geometrically noisy point cloud. 






Fig. 7.11. Comparison of our meshing scheme and the Tight Cocone algorithm EF on a ge¬ 
ometrically noisy point cloud. All points are considered in the computations. Left: Our meshing 
result with a zoom-in of the nose. Right: The result of the Tight Cocone algorithm with a zoom-in 
of the nose. 


for a quasi-conformal map in [6] for achieving better conformality. Thirdly, we intro¬ 
duced a balancing scheme for guaranteeing an even distribution of the spherical point 
cloud parameterization. Experimental results show that our proposed algorithm is 
highly efficient and accurate. With the aid of the spherical conformal parameteriza- 
tions, almost-Delaunay triangulations and high quality quadrangulations of genus-0 
point clouds can be effectively created. The meshes generated are guaranteed to be 
of genus-0 and no post-processing is needed. Besides, our meshing method is sta¬ 
ble under geometrical and topological noises on point clouds. Moreover, multilevel 
representations of genus-0 point clouds can be easily computed. As a remark, our 
proposed spherical conformal parameterization algorithm also works efficiently on tri¬ 
angular meshes. In the future, we plan to establish a rigorous theoretical proof of 
the convergence of our parameterization scheme, and extend our method to handle 
disk-type point clouds and point clouds with arbitrary topology. 
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Fig. 7.15. Multilevel representations of a genus-0 lion vase point eloud. For better visual¬ 
izations, we create mesh structures on the representations. Top left: a lion vase point cloud with 
256094 points. Top middle to Bottom right: the multilevel representations with 0, 1, 2, 3 and 4 
subdivisions. The representations are with 696, 2778, 11106, 444^^ 177666 points respectively. 
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